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Abstract. We prove global existence and uniqueness of solutions of Oldroyd-B systems with 
relatively small data in K , , in a large functional setting (C a D L 1 ). This is a stability result, 
solutions select an equilibrium and converge exponentially to it. Large spatial derivatives of the 
initial density and stress are allowed, provided the L°° norm of the density and stress are small 
enough. We prove global regularity for large data for a model in which the potential responds 
to high rates of strain in the fluid. We also prove global existence for a class of large data for a 
didactic scalar model which attempts to capture in the simplest way the essence of the dissipativc 
nature of the coupling to fluid. This latter model has an unexpected cone invariance in function 
space that is crucial for the global existence. 



1. Introduction 

The complex fluids models we consider treat the interaction of a normal fluid (incompressible in 
our case) with insoluble matter suspended in it. Models have been devised to deal with microscopic 
elastic thread-like objects such as polymers ([12]. pfj]b The complicated hydrodynamic interactions 
are simplified using the separation of scales, replacing the many degrees of freedom due to them 
by few representative ones. In the models we consider here, an end-to-end vector in R 3 (m in this 
paper) represents the orientation of the polymer and it is the sole variable retained to describe the 
microscopic object. The models that are most studied with this degree of simplification include 
the kinetic description of the Oldroyd B model, the FENE model and variants. The probability 
distribution of the vectors m depends on time and physical space, and it is assumed to be absolutely 
continuous with respect to the usual Lebesgue measure dm, so the chance of finding m at time 
t and location x € R 3 in the volume element dm is f(x,m,t)dm. The aim is to describe the 
evolution of / in space and time. When the polymer concentration is so dilute that the polymers 
do not interact, and when the fluid is smooth, then the equation of / is a linear Fokker-Planck 
equation. The effect of the fluid can be understood perturbatively: because the objects suspended 
are assumed to have faster time scales and shorter length scales than the scales on which solvent 
varies, it is then justified to treat the fluid as varying little as far as the microscopic suspended 
objects are concerned. This macro-micro interaction can be rationally discussed, and leads to a 
kinetic description of the particles, given a fluid flow. On the other hand, the collective effect 
that the particles have on the flow itself is a macroscopic effect of microscopic insertions: a micro- 
macro interaction. These interactions are much more mysterious. There have been attempts to 
produce systematic upscalings based on non-equilibrium thermodynamics formalisms ([2], |13j). 
When the microscopic insertions have a larger number of degrees of freedom, m € M , where M is 
a Riemannian manifold representing a finite number of degrees of freedom with constrains, then the 
kinetic equation is a Smoluchowski equation on the manifold. The upscaling principle advocated 
by one of us ([5], [5]) is easiest formulated as the requirement that the sum of the energy of the fluid 
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and the free energy of the microscopic quantities be a Lyapunov functional for the coupled system. 
In the known examples, this requirement leads to familiar rules of determining the added polymeric 
stress from the micro- micro and the macro-micro interactions (|16J). The mathematical treatment 
of the coupled systems is far from being complete but has received recently more attention. Early 
work ([2], [IE]) established local existence results for Oldroyd-B and FENE type equations. The 
problem of existence of weak solutions is rather open. Global existence of weak solutions in the 
presence of spatial diffusion of the polymers was proved in a sequence of papers, the most recent of 
which is ([I]). Global existence of weak solutions via propagation of compactness was proved under 
the corotational assumption ([22], [23]) and, very recently, for the full FENE model ([25]). There 
is no such result for the Oldroyd B model. The global existence of smooth solutions for small data 
for Oldroyd B-type models was established in ([T7], [S])- Global existence of smooth solutions 
for large data in 2D was established for Smoluchowski equations on compact manifolds ([5J, [7J, 
[3], [TU], [H], [37J ). Global regularity for large data in the FENE case, under the corotational 
assumption was proved in Ql9|. |24j ) . An approach based on Lagrangian particle dynamics was 
described in ([20]). Sufficient conditions for regularity in terms of bounds on the added stress tensor 
were established in ([3], [T5]) and further refined in ([H]). Numerical evidence for singularities was 
provided in ( |30j ) . 

In this paper we address issues related to global existence of smooth solutions in simplest kinetic 
models. In the presence of a quiescent solvent, the polymer distribution is the unique stable time 
independent solution of the linear Fokker-Planck equation 

d t f = ediv m (/V m (log/ + [/)) 

i.e. / = Z~ 1 e~ u with Z a normalizing constant. The simplest form of U is harmonic, U(m) = ^jj- 
and then the solution is Gaussian. The constant R 2 represents the expected value < \m\ 2 > of the 
square of the length of the end-to-end vector. In the presence of smooth fluids the kinetic equation 
changes and acquires a dependence on the macroscopic independent variable x. 

In the second section of the paper we provide a priori bounds for linear Fokker-Planck equations 
with Oldroyd B type potentials. We do this by essentially giving a formula for the solution of the 
Fokker-Planck equation in terms of the particle trajectories of the underlying flow. 

In the third section we derive a classical estimate for the velocity gradient in terms of the added 
polymeric stress, when the coupling to fluid is done via time-independent Stokes equations. 

In the fourth section we prove global existence results for small data in C a , at arbitrary Deborah 
numbers. The previously known global existence results for small data are done using energy 
estimates and require derivatives of the stress. The main difficulty in obtaining bounds for small 
data is the fact that quantities are not monotonically decreasing, and, in particular, the spatial 
gradient of the density of particles can grow in time, but it nevertheless saturates. The system 
formed by the kinetic equation for the particle distribution and the velocity of the fluid can be 
reduced to a closed system for a reduced added stress r and the particle density. This system 
looks like a damped and driven Euler equation where r (a matrix) plays the role of vorticity and 
/?, the transported density is part of the driver. The global existence of solutions is proved under 
the assumption that the L°° norm of the initial data is small enough (| 1 T3[) . The class of initial 
data that lead to global existence includes densities and stresses with arbitrarily large gradients, 
provided the L°° norms are suitably small. While the limit added stress r and limit velocity vanish, 
the limit particle density is a re-arranged initial particle density, a deformation of the initial density 
along the limit back-to-labels map. The class of solutions we discuss is in fact rather wide, and 
in particular the gradients need to be understood in distribution sense. The proof of uniqueness 
of solutions is done using Lagrangian transformations. An Eulerian proof is not readily available 
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at this moment. In order to explain this we pursue the analogy with Euler equations. Both the 
damped Euler vorticity equation and the equation for r have the form 

D t r = -ct + N(t) 

where c > is the damping and the nonlinearity N is quadratic and non-local, obtained from 
products of t and principal value singular integrals Mr. Here D t is material derivative along the 
divergence-free velocity u. The local (and global for small data) existence results are in Banach 
algebras in which TV is a continuously differentiable nonlinear function. One of the largest and most 
natural such Banach algebras is C a , < a < 1 . If we consider the difference S of two solutions t\ , 
t 2 , it obeys the equation 

D t S = -cS + N'(t)S - v ■ V x t 

where D t is transport along the average velocity u = \ {u\ + u 2 ), f is the average of the solutions, 
t = i(ri + t 2 ), N'(t) is the linearization of N obtained using polarization, and v is the velocity 
difference v = u\ — u 2 - The last term in the equation is present because the velocities depend 
on t: mi is obtained from n and u 2 from r 2 . In both the (damped-driven) Euler case and in the 
Stokes-Fokkcr-Planck system we study, the dependence of u on r is linear. The term v depends 
therefore linearly on 5. Uniqueness proofs based on energy estimates proceed by estimating the 
growth of 5 in some norm. As the initial data vanishes, any closed, super-linear estimate of the 
norm is sufficient to prove uniqueness. (We could allow logarithmically sublinear estimates, but 
that is as far as that can go, in general). There are two sources of difficulty in obtaining these 
estimates. The transport term u ■ V X S can be dealt by using integration by parts, if the norm of 
S we consider is an L p norm. The other source of difficulty is V x t. If we are in the framework of 
C a n L 1 adopted in this paper, then v is very nice but in order to deal with V x f we need to work in 
H~ x . In general, the presence of both v ■ V x t and u- V x <5 together is lethal for an Eulcrian, energy 
method approach, unless we work in a space in which V x t is tame, like for instance IU 1 ' P (R <1 ), with 
p > d. The reason there is a successful Eulerian proof of uniqueness for Euler equations for the 
vorticity in C a is that, for Euler equations, the vorticity equation comes from a velocity equation 
with a good cancellation property. Thus we can "retreat" to one less derivative in the equation for 
S. In terms of the vorticity equation this means integrating against a vectorial stream function, 
i.e. working in H~ l . In the Euler equtions the transport term together with the term N'(t)S 
can be integrated by parts one more time, revealing only H~ l terms in <5. In the absence of such 
algebraic reasons, the retreat to one less derivative does not appear to work, and the fix for one 
difficult term does not work for the other, and vice-versa. The Lagrangian approach introduces 
commutators that are well behaved, and is successful. 

In the fifth section we show that the system is regularized for large data if the potential is 
allowed to respond to excessively high rates of strain S in the fluid. The Oldroyd-B potential is 

i2 

harmonic, U = as it were generated by a spring force that is related to a material property 
of the polymers. But, in fact, the representation of the polymer by m is over-simplified, and R 
has something to do with an average restoring force in the ensemble, and not with the maximum 
extension allowed for an individual molecule. In the case there is feedback coupling to a fluid, the 
interaction might have properties that depend of the local properties of the fluid. We show that, 
if we allow R to depend on the local rate of strain in the fluid, when this exceeds a threshold, 
and to grow with it, = S(S), then the equation has global smooth solutions for all large data. 
Here 5(S) can vanish if the rate of strain S is not too high, and becomes asymptotically linear in 
S for high values of (norms of) 5. Allowing R to grow in response to high rate of strain is like 
allowing more entropic "slack" in the molecules at high rates of strain, and this turns out to be a 
mathematical regularizing mechanism. 
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In the sixth section we return to the Oldroyd B large data problem, and give an example of blow 
up in one dimension for a crude model that does not respect the principle of free energy decay. 
We motivate and give then a simple scalar example of an equation that has some of the features of 
the type of equations arising in the large data Oldroyd B equation, and which has an unbounded 
set of initial data that lead to global smooth behavior. The reason for the global existence is an 
invariance of a cone in function space under the nonlinear evolution. The consequnece of this 
invariance is an a priori bound that is sufficient for the persistence of smoothness of solutions. 

2. Linear Fokker-Planck Equation: kinematic observations 

We take a vector field u(x,t) and a scalar f(x,m,t) representing a two or three-dimensional 
incompressible velocity and a particle distribution. We start by describing the simplest particle 
distribution equation 

(1) D t f + {{V x u)m) ■ V m f = ediv ro (/V m (log/ + U{m))) 
Here D t — dt + u ■ V x is material derivative. The potential is given by 

(2) U im) = feg. 

In this section R is a positive constant. We consider here the case m £ M d . We associate to the 
particle distribution / an added stress tensor: 



(3) a(x,t) = / (m®V m ?7) f(x,m,t)dm 

JmeK d 

In this normalization R has units of length, e/R 2 is an inverse time. The potential U and the 
stress a are nondimensional. We study first (JTJ), without regard to the coupling to the fluid. In 
what follows, the velocity field could be quite arbitrary, but it is assumed to be smooth enough 
for the calculations below: the variable x £ M. d , and the velocity field u(x,t) is in C 1,a (K ), is 
divergence- free 

div^u = 

and decays at infinity. Before we start our computations, let us make a few general comments 
regarding the matrix a. First, because the potential U is radially symmetric © and from its 
definition it follows that a is symmetric and non- negative. Moreover its off-diagonal entries 
are bounded by the trace 

(4) \aV(x,t)\<±Tr(*(x,t)). 
We take now the Fourier transform in m 

(5) f(x,t,t)= I e- m <f(x,m,t)dm 
In view of fl2J, / obeys 

(6) D t f(x, £, t) + [JLl - (V x uf] £ • V e / = -e\e?(x, £, t) 

We used the fact that u is divergence-free and denoted (V x w) T the transposed matrix. We solve ^ 
on characteristics. The connection between the Lagrangian paths of u and the Oldroyd B equation 
was exploited in [20] in the less complicated case e = 0, p = 1. Let X(a,t) be the particle paths, 

(7) d t X(a,t) =u(X(a,t),t) 
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with X(a, 0) = a and let 

(8) g{a,t) = (V x u)[X(a,t),t) 
Let now £ = £(a, rj, t) solve the ODE 

(9) ^=^. g{a , tn 

with initial data £(a, 77, 0) = 77. We take the fundamental matrix $(a, t), solution of the linear ODE 
system 

(10) ^(a,t) = -g(a,t) T $(a,t) 
with initial data <£>(a, 0) = I and we have then 

(11) Z(a,r ] ,t) = elF$(a,t)r ] 
We write now 

(12) F(a,r,,t) = f(X(a,t),Z(a,ri,t),t). 
Then ((6|) implies 

and integrating we obtain 

(13) F(a, ry, t) = e- £ /«.' I^^)! 2 ^/^ ) 
Now we invert the linear map 77 H ► £(a, 77, t), and write from (jlip. 

(14) » ? (a,e > t) = e-^*(o,t)^ 
where ^(a, t) = $(a, obeys 

(15) ^*(o,t) = *(o,t) ff (o,t) T . 

with initial data ^(a, 0) = I. Reading (IT21 at 77 = 77(0, for a fixed £ and substituting in ([T^]) 
we obtain 

(16) f(X(a,tU,t) =e-^o\6^Ma,^)\ 2 ds^ e -^ MO 

where fo(x,£) is the Fourier transform in m of the initial data fo{x,m) = f(x, m, 0). Let us 
consider now the "back-to-labels" map A(x,t), inverse of X(a,t). Let 

(17) M(x,t) = V(A(x,t),t) 
and note that it obeys the transport equation 

(18) D t M = M{V x uf 

with initial data M(x, 0) = I. Reading JT6]) at a = A(x, t) and we deduce 

(19) f(x, £, t) = e~ e ti \e(Mx,t)rtA(x,tU,t),s)\ 2 dsf o ( A ( x ^ t ) je -$M[x, t)C). 
We compute directly from (fTTj) and (fl4|) 

(20) £{A(x,t),Ti(A(x,t),S,t),s) =e-^Q{x,t,s)t, 
with 

(21) Q(x,t,s) = q(A(x,t),t,s), for /j > s, 
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(22) q(a,t,s) = $(a,s)^(a,t) for t > s, 

and deduce thus from (fH 



(23) f(x,U) = e- e SS'~^^W^')^%(A(x,t),e-^M(x,t)0 

Let us note now that the matrix a 13 (x, t) is computed from the Hessian of / at £ = 0, 

(24) ^) = "i^^%=° 
Using (|23p and noting that the cross-terms vanish, we deduce 

= ifpOM) Joe" 2 ^ [Q T (a:,t,a)Q(a:,t,a)] ij .da 



2rt 



-^-M kl (x, t)M l3 (x, t) £4^ (A(x, t) , 1 ?=0 



We denoted by p(x,t) = f(x,0,t) — L d f(x,m,t)dm. In view of JT]) this obeys the transport 
equation 

(26) D t p = 
and therefore it is given by 

(27) p(x,t) = f (A(x,t)) = p (A(x,t)) 

in terms of the initial particle density at x. We used this in (|2~5j) as well as the summation 
convention. In terms of the initial stress, the expression is 

ff«(x,t) = %p{x,t)^e- 2 -^ [Q T (x,t,s)Q(x,t, s )] i3 ds 
+e~WM kl (x, t)Mij(x, t)a kl (A(x, t),0). 
The solution in ([28)) solves the equation 

(29) D t <r = -|jff + (V x u)a + cr(\7 x u) T + ^p{%, t)l 

which can be easily derived from the equation (fTJ by multiplying with jpmirrij and integrating 
dm. Let 9(x,t) be a passive scalar, i.e., a solution of 

(30) D t 6 = 0. 

Then, if M solves (|18|) then MV x 8 is again a passive scalar, i.e., 

(31) D t (MV X 9) = 0. 
This can be easily checked because 

Dt{dj0) = -(d jUl )di6 

and so 

lh<M,,i)J)- = M lk {d kUj )d 3 e - M l3 id jU i)di6 = 

This means that if 6(x,t) = do(A(x,t)) with arbitrary smooth 9o(a) then M(x,t)W x 8(x,t) = 
(V a 6*o)(A(x, t)) holds. In fact, because the initial datum of is the identity matrix, and because 
of uniqueness of ODEs, it follows from (|15|) that 

(32) *(o,t) = {V a X) T (a,t) 
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and consequently 

dX i 

(33) M ki (x,t) = — (A(x,t),t). 



From (J28J) and ([33]) it follows that that 

/ 34 n (r(z, t) = j?p(x, t) e l^Q T (x, t, s)Q(x, t, s)ds 

+e-^(W a X(A(x,t),t))a Q (A(x,t))(W a X(A(x,t),t)) T 

Let us introduce the notations 

(35) e\x,t) := (WX l )(A(x,t),t) 

(36) e l3 { Xl t) := <To(A(x,t)y(x,t) ■ e J (x,t) 
Note that 

(37) a{x,t) = ^p(x,t) J e- 2 -^Q T (x,t, S )Q(x,t, S )ds + e~l$ (e ij (x,t)).. 
Moreover, 

(38) D t e* = (V x u)e l 

with initial data that are constant in space and equal the canonical basis of R d , e l (x, 0) = (Sij)j 
and that q given in ([22)) is given in terms of the gradient V ' a X by 

(39) q(a, t, s) = [(V„X(a, s))~ l ] J '(V ' a X{a, t)f 
solving 

(40) d s q(a 7 1, s) = -g T (a, s)q(a, t, s) 
and 

(41) q(a,t,t)=l, q(a,t 7 0) = (V a X(a,t)f. 
Passing to Lagrangian variables in (|37l) we obtain 

u 2 ) a(X(a, t), t) = jpp {a) / * e ^g T (a, t, s)g(a, i, s)ds 

+e-^ (V X(a, t))tr (a)(V a A-(o, i)) T 
and integrating by parts we deduce 

(43) CT ( X M),*)-£oJ>)II = 

2 Pa (a) /„' e -TV(a, t, s)S(X(a, s), «)?(o, t, s)ds + e"^ (d a X(a, *))r (a)(d a X(a, t)) T 
where we introduced the reduced stress 

(44) r(x,t) = a(x,t) - p(x,t)l 
and where 

(45) S(x, t) = i [(V x u(a;, t)) + (V^^, t)) T ] 

is the rate of strain. Returning to Eulerian variables, we have 
r(x,t) = 

(46) 2p(x,t) f* e'^^ Q T (x,t, s)S(X(A(x,t), s), s)Q{x,t, s)ds 
+e-W (V a X(A(x, t),t))T (A(x, t)){V a X(A(x, t), t)) T 
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It is clear that in general, once the velocity u(x,t) is given, we can compute everything we need 
to know about r (and hence a) from (|46[) . We will be interested in particular in the relationship 
between various norms of u and norms of r, as these relationships will serve in establishing bounds 
for solutions of nonlinear equations. It is clear that we should start with the fields 6j, and in order 
to understand their norms we perform some standard Lagrangian estimates. We take the finite 
difference 

(47) 5 h (d a X\a, t)) = d a X\a + h, t) - d a X l (a, t) 
Its equation follows from ([7]) by differentiation, 

(48) d t (W a X l ) = g(a,t){V a X i ) 

and then taking the finite difference: 
(49) 

—d h [daX (a,t)) = 6 h [d a X (a,t)) + dhg(a,t) 

where g(a,t) is the velocity gradient in Lagrangian coordinates (JSJand 
(50) S h g(a, t) = g(a + h, t) - g(a, t). 



For matrices L we use the notation \L\ for the Euclidean norm of the matrix \L\ — \/TrL*L. We 
denote 

(51) 7(t) = sup|g(a,s)| = \\V x u(-,t)\\ L ^ {dx) 

a 

a quantity of some importance in the sequel. From P5|) and the fact that the initial data for d a X l 
is {5ij)j=i....d we have that 

(52) sup|<9 a X>,i)| < e ^^ s)ds 

a 

Because the initial data for 5h [d a X % ^ vanishes, we have, from Gronwall's inequality and (1521) : 

(53) \S h (d a X*(a,t))\ < e/o-rM* / \S h g(a,s)\ds 



We are interested in quantities in their Eulerian form. We consider the Holder seminorm 

54) [4> a = sup — : — 

x ^ v \x-y\ a 

with < a < 1. We note that the back-to-labels maps are Lipschitz and 

(55) sup M*,t) ~ A{y,t)\ £eJoS( , )d , 

because 

(56) D t A = 
with initial data A(x, 0) = x, and consequently 

(57) D t V,A = -(V x A)(V x it). 

In view of (|55|) and the fact that Holder seminomas behave nicely with respect to compositions 



with Lipschitz functions, 

[4>oA] a < [4>] a \ a 
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if A is a bound on the Lipschitz seminorm of A, it follows that 

(58) [e l (;t)] a < [VI'(.,()]/^( s ) ds 
and, similarly 

(59) [g(;t)] a < [V x u(;t)] a e"ti^ ds 
From these considerations and from (1531) it follows that 



(60) [e'M)] a <e( 1+2a »£^ ds [V x u(;s)] a ds 



The reason for the 2a "loss" in comparison to the Lagrangian estimate (|5"3")) is that we "loose" twice 
by composing with Lipschitz functions, once because the objective of the estimate is composed, 
and another, because what what we estimate it with is also composed. We also obtain, in the same 
manner 

(61) \\V x e l (;t)\\ LP{Rd) < e 3 /o-rW^ f \\V x V x u(;s)\\ LP(m ds 

Jo 

We perform similar calculations for the matrix Q(x, t, s) taking into account (|2T1 l22l HQl ST}. We 
deduce 

(62) ||Q(-,M)||l~ ( r-) <Ce/.'^ 



(63) \S h q(a, t,s)\< Ce?. ^ z)dz \ \S h g(a, z)\dz 

J s 

and consequently 

(64) [GO, *,*)]« < Ce^ 1+2a ^o^)^ f [V x u(-,z)] a dz. 



with C = y/d depending only on the dimension of space d. We are ready to state the estimates on 
t in terms of u: 

Proposition 1. Let u be a divergence-free function belonging to L 1 (0,T;C 1 ' a (M. d ). Let o~o(a) be a 
Holder continuous, L 1 (M. d ), positive symmetric matrix and let po be a positive, Holder continous, 
L 1 (M. d ) function. Then t = a — pi, given in the expression j46\ ), obeys 

IK-,*)llx/»(R<*) < C\\po\\l°°(bl*) fa expl-^r^ + 2 j(z)dz\j(s)ds 
(65) r , S J 

+C\\r \\ L oo m exp [-^ +2f* ~f(s)dsj, 

and 

H;t)] a < C[po] a e a fo^ dz J t exp{-^^+2j t sl (z)dz} 1 (s)d 
+C\\p \\ L ~ {Md) /* exp \- 2 -^ + (2 + 2a) f* j^dz}^) f s [V x u{; z)} 



(66) 



s 

dzds 

C||Po||l°°(r") Jo exp |- 2e( ^ s) +2 f s j(z)dz + 2a J* j(z)dz\[V x u(-, s)] a ds 
[ro]a + \\to\\l°°(r*) fo[Vxu(',s)] a dt exp{-fp + (2 + 2a) /„* 7(a)ds| 



-C 



with 7(f) = sup^, |Va;u(x, f)| and the constant C > depending only on the dimension of space d. 
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Proof. The L°° estimate follows from (|52"j) and 

(67) \W{-,t)\\L^ (dx] <e^^ ds 

which follows immediately from (155)) and the fact that the initial data are of unit norm. The Holder 
seminorm estimate follows from the algebra inequality 

[#]a < MalMU- + II^IU-Ma: 

(|62|) and (l67|) above, (|60|) , (|64|) , and the behavior of the Holder continuous function under com- 
position with A and (|55|) . The VF 1,P estimate follows using the product rule. This concludes the 
proof of the proposition. 

Before we conclude this section, we comment on the case of d = 2. In this case the velocity 
gradient matrix is 

(68) V.«(M)=( , fi%, t) ^JT) 

yn{x,t) + - ! ^ Lj - -X(x,t) j 

where 

(69) X(x, t) = - (d 1 u 1 (x, t) - d 2 u 2 {x, t)) , 

(70) n{x, t) = ^ (9iu 2 (s, t) + d 2 v} (x, t)) 
and 

(71) lu(x, t) = d!U 2 (x, t) - d 2 u 1 (x, t). 

The symmetric part of the gradient (rate of strain matrix) is 

X(x,t) H>{x,t) 
n(x,t) —X(x,t) 



(72) S(x,t)-- 
Differentiating ([55)1 we can write 

(73) D 2 e l = ({V x u) 2 + D t (V x u)) e l 

This equation holds in any number of dimensions, but because of the Hamilton- Cay ley theorem, 
in two dimensions we have 

(74) (V x u{x,t)) 2 =S(x,t)I 
with 

(75) S(x,t) = X 2 (x,t) +fi 2 (x,t) - ^uj 2 (x,t) = -det(W x u)(x,t). 
Therefore the equations (f?3"j) become 

(76) D\e l = (SI + D t (V x u))e l 
Similarly, we have 

(77) D 2 s Q = -Q(SI+D s (V x u)) 

We can easily integrate the equations (|76[ [77| when the velocity gradient is constant in space and 
time (so that D t (V x u) = 0). In this case 6 is a constant and its sign dictates the size of solutions 
e'and Q. Consequently, when S < or when S > but y/6 < then r is bounded, and when 
5 > (^) 2 j then r may grow exponentially in time. Integrating (l76l [77l) when D t (V x u) = we see 
that, if S < the functions e l and Q are bounded in time and when 5 > 0, they grow like exp \/5t 
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in time. If the initial distribution fo is radially symmetric in m then the matrix ao(a) is a multiple 
of the identity matrix (a scalar matrix, or an isotropic stress), 

a (A(x,t)) = cp (A(x,t))l. 

If 8 is a positive constant, integrating (|T6]) we have 

i / cosh(tVI) \ 2 ( f) _f sinh(tVS) 

6 [X ' t} ~ \ sinh(t^) ) B [X,V ~ { co S h(tVS) 

In this case _ 

e u (a;,t) = e 22 (x,t) = cp Q {A(x, t)) cosh(2tVS) 

e 12 {x,t) = cp {A(x,t))sinh(2ty/S) 
The case of a uniform gradient is used in rheology, but if the velocity gradient is square-integrable 
in space, then 5 has to have average zero in space because 

(78) (V x n)?- = d k (u i d j u k ). 
Therefore, for moderately varying Vi« we can expect r to be bounded. 

3. Coupling to the Stokes system 

The Stokes system is 

(79) — A x u + V x p — kdiv x a 
coupled with 

(80) div^u = 

The coefficient k has units of inverse time, and represents the ratio between an energy (kT) per 
unit mass and kinematic viscosity. The fluid density is taken to be one. We discuss the situation 
in which the fluid occupies all M. d and is at rest at infinity. The velocity and pressure decay in 
space and are obtained then from a by classical singular integrals. We can equally consider the 
case in which the velocity and pressure are periodic in space. Let us denote by 

(81) Rj = A- 1 ^- 
for j = 1, ... d, in R d the Riesz operators where 

(82) A=(-A)2 

is the Zygmund operator. By modifying the pressure, we may write 

(83) - A x u + V x p = kdiv x (T) 
and using (|8"D]) to solve for p 

(84) p = -kR m R n (r mn ) 
we deduce 

(85) u l = kk- 1 (^(r li ) + R i R m R n (T mn )) 
and 

(86) dju* = kRj (R(t u ) + R i R m R n (r mn )) 
which we abbreviate as 

(87) V x u = kllr 
The system formed by (|3"9l H3")l together with (|55|) is closed. 
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3.1. Energetics. The system formed by (p} coupled to (1791 has a Liapunov functional. We take 
(HJ), multiply by fc(log/ + U), integrate, take (JTSJ) multiply by u, integrate, and add the two: we 
deduce 

(88) k Tt Ur' /*■ / lo § f dmdx + 3 /*«■ ^(as, *)]<**] 

= — Jr^ iV^wt^, i)| 2 da; — fee J Rd / Rd / |V m (log/ + Z7)| 2 d/TzcZx 

This is a general property of the class of Fokker-Planck equations we are interested in, but in the 
case of the Oldroyd B system we are considering here, we can obtain the energetics directly from 
the equation ([29]) by taking the trace, adding to the Stokes equation (|79l) multiplied by |w and 
integrating: 

(89) d t [ Tradx + \( \\7 x u\ 2 dx = f Tr adx + || f p dx 

JR d & JR d n JR d JR d 

Using the fact that Tra together with p bound the entries in r, we derive 

(90) ||r(-,t)|| L i( R d) < C[||r ||ii(Hi) + ||A)||i 1 (E' 1 )] 
with a constant C depending only on d. 

3.2. Gradient bounds. We are going to investigate bounds on the time integral of the maximum 
gradient of velocity, 

rT 

\\V X U(-,t)\\ L oa (dx) dt 

In order to do so we recall a classical fact ([29], Chapter 3, section 5), namely that operators given 
at Fourier level by multiplication by functions that are homogeneous of degree zero and are smooth 
on the unit sphere can be represented as sums of multiples of the identity and singular integral 
operators of a classical type. 

Proposition 2. There exists a constant C depending only on d and a such that, for any a, 

a d 



(91) \\na\\ L ~ (dx) <C\\a\\ L ~ (dx) { \ - ]«,. 

where 



1 



\\0-\\L°°{dx) 

(x)-a mn (y)\ 



(92) [a] a = maxsup- 

mn x^y \x - y\ a 

and < a < 1 . There exists a constant C depending only on d and a such that 

(93) [Ka] a < C[a] a 

Proof. In view of the fact that TZ is made from operators of the form RjRk and RjRkRi we have 
the representation 

{K3)ij{x) = C ijmn a mn (x) + P.V. f h ^ y) a mn (x - y)dy 

jR d \y\ 

where Cij mn are constants and each kij mn is a smooth function (actually a harmonic polynomial 
of degree less or equal to three) of y = y/\y\, with mean zero on the unit sphere § d_1 . We need 
to prove the inequality for each entry, and we reduce the proof to the scalar case. The proof then 
follows along very classical lines. We break the integral 



Ka(x) = P.V. [ ^a(x - y)dy 

jR d \y\ 
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into three pieces. We choose two numbers m < M and write 

h{x) = P.V. 



Kv)~i 

y\<m \V\ 



and 



h{x) = 

h{x) = 



m<\y\<M W 



■jcr{x - y)dy 



— ja(x - y)dy. 



'\y\>M \V 

We have, with the cancellation property of k(y), 

\h(x)\ < Cm a [a] a 

then clearly 

M 



\h(x)\ <C\og — ||of|| L o 
m 



and 



\h{x)\ < CM- d \\a\\ L i (dx) 



Now we choose the two length scales, 



and 



M 



<r\\LH 



\a\\L°°{dx) 







2 + log 


(-)} 




\m j 



If M > m we use both in the inequalities above and we obtain 

\Ko(x)\ < C\\a\\ L ^ (dx) 

If M < m then we forgo the division above and split the integrals into only two pieces, up to M 
and from M. We obtain then 

\Ka(x)\ < 2C\\a\\ L ^ dx) 

This ends the proof of (|91l) . The inequality is a classical inequality for singular integrals of 
the type above. This ends the proof of the proposition. 

4. Small data 

The approach we choose to control [r] a is via the explicit formula (|46|) . The equation obeyed 
by r can be derived directly from 

(94) 



D t r = -|jr + {V x u)t + t(S7 x u) t + 2pS 



Because 2S = klZr + fc(7?.r) T , depends linearly on t, it is clear that strong enough damping will 
prevail, sending the solution to zero. It is clear however that 



(95) 



> C||p |U°°(R<*) 
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will have to be part of the requirement for this stability result. The system of ((T|) with @ and ^ 
coupled with ([75)1 is equivalent to the system formed by the Stokes system (|83p and (f2lj)). We 
use Proposition ([T]). We first prove a good local existence result. Let us denote 

(96) Mi = ||po||ii(Rd) + UtoIIli^) 

(97) Moo = ||polU~(Rd) + ||To|!Loo( R£ i) 
and 

(98) M a = [p ] a + [t ] q 
Let us introduce the Deborah number 

(99) 0:=- = —, 

k e 

Ft 2 

the nondimensional ratio of the particle time scale — to the viscous fluid's response time to the 
stress added by the particles, We denoted half the damping rate by 

(100) Ko - jp- 

Theorem 1. Let (r ,po) € L 1 (R d ) 2 n C a (R d ) 2 . There exist constants e > 0, V > 2 and there 
exists a time T > and a weak solution of pffi $5) (p,r) £ C([0, T ], W~ 1 ' p (M d ) 2 ) n 
L°°([0, To]; (C a (R d ) n T 1 (R d )) 2 ) ; 1 < p < oo, satisfying the equations in weak sense, and such 
that 

(101) DkoToMoo |l + log (l + M^M^M~A \>e 
and 

(102) \\p(;t)\\^m + M;t)\\Li(R') < TMx, 



(103) ||p(-,«)IU-(r«i) + lk(-,0IU-(R*) < rM °° 

and 

(104) W',t)]a + M-,t)] a <M a 
hold on [0,T ] 

Proof. We consider the system (|94[ 1551 |21)|) with initial data that have been regularized by 
convolution with a standard mollifier <pg. We obtain uniform bounds for quantities of interest and 
pass to the limit 5 — > 0, removing the mollifier. In order to simplify the exposition, we will denote 
the solutions t^ s \ p^ by r, p as if they did not depend on 8. The bounds and constants below will 
indeed be (5-independent. We use the following notations 

(105) 7(t) := {Dko)- 1 sup || V x u(-, a)|| L » (Rd) , 

0<s<t 

(106) 6(t) := (Z?^)- 1 sup [V^(-,s)] Q , 

0<s<t 

(107) n(i) := sup [||r(-, s)||ioo (E <j) + \\p(-, s)!^^)], 

0<s<t 
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and 

(108) m(t) := sup {[t{; s)] a + [p(-, s)] a } 

0<s<t 

We prove first uniform bounds, then discuss their use. We use ([90)) 

(109) sup [||r(-,s)|| L1(Rti) + \\p(;s)\\ L i m < CM 1 

0<s<t 

which is based on the fundamental dissipative property of the system. We note that the function 
(x, a, b) \-t x log (l + 22) is increasing in each of its arguments x > 0, a > 0, b > 0. Therefore, we 
may use upper bounds in the right-hand side of (|91l) . Using (I91[) we deduce 

(HO) 7W<an( t ){l + log( 1 + ""»y )} 

and using 

(111) b(t) < Cm(t). 
We proceed to estimate in (p5|) : 

(112) n{t) < CAfoo[l + Dn tl(t)}e CDKot ~^ 
and in (|66[) 

(113) m(t) < CM a [l + DK a tj(t)]e CDKot ^ t) + C(M oo £>K t)[l + DK Q tj(t)}e CDKot ^ t} m{t) 
The term DKotj(t) can be bound from (|110[) 

(114) D Ko tj(t) < (CDK t)n(t) |l + log (l + ^^^ ^j | 

The structure of the inequalities is the following: if we denote the group 

(115) x(t) := DK Py(t), 
we have from (|112[) 

(116) nit) < CM^l + x{t))e Cx ^ 
and therefore, from (|114[) we deduce 

x(t) < (CDkoW^Q. + x(t))e c ^(l + log(l + cil+x{ l ))e a^ )) 
h + log (l + m{t)^M^M-A } 



(117) 



x 



so 



(118) F{x(t)) < C{M 00 DnQt) {l + log (l + m(t)& Mf+ a M' 1 ) } 
where 

(119) F{x) = xe- c *(l + x)-i (l + log(l + c{l + x)eCx 
Let 

(120) a{t) := MaoDnat 
Now (|113p can be written as 

(121) m{t) < C(l + (x(t))e Cx{t) [M a + a{t)m(t)} 
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and 
(122) 



F(x) < Ca(t) {l + log (l + m(t)^ Mf+* M~ l ) } 



We fix Co a large enough absolute constant, larger than the C we encountered so far. 
loss of generality we may assume 1 + x < e Cox and 



Without 



1 



C(l + x)e 



so that 



F(x) > xe 



-2C x 



c a 



and therefore, as long as x(t) < 1 we have from (|122j) 



(123) 
and, from (fT2Tj) 

(124) 

Using x(t) < 1 and requiring 
(125) 



x(t) < C 2 e 2Co a(t) {l + log (l + m{t)T^M*+ a M' 1 ) } 



m(t) 

m7 



1 + a{t) 



m(t) 



a(t) < iCJ-V** 



we deduce from (I124|) that, as long as < 1, it follows that 

m(t) 



(126) 



< 2C e 



2C 



Returning to a;(t) we have that, as long as x(t) < 1, it follows from f|123[) and (|126l) that 

(127) x(t) < Cle 2C °a{t) |l + log (l + (2C a e 2C ° M a )^ M? 
We require therefore also 

(128) a(t)<±{c 2 e 2C °{: 



1 + log ( 1 + (2C e 2Co M Q )^M 1 d+Q M, 



'»)} 



and thus, as long as < 1, it follows that actually x is bounded by x(t) < i. The initial datum 
for x(t) is zero, and a(t) is explicitly proportional to t and starts from zero. Therefore, because of 
continuity, if the time t is taken small enough so that a(t) satisfies the requirements (11251 1128|) . it 
follows that x(t) < j on that interval, because we can reason by contradiction, and the first time 
it would reach | it would have to be not larger that |, which is absurd. Consequently, it follows 
from ([T2"B]) that m(t) < 2C e 2C °M a and from (fTig|) that n(t) < 2C e 2C °M oo on the time interval, 
providing the required short time uniform bound. The rest of the proof of the theorem is based on 
the uniform bound. It is well-known that regular solutions exist as long as the L°° bounds persist 
(see [3], [IB], also next section). Therefeore classical solutions with regularized initial data exist 
for this interval of time. Indeed, by regularizing the initial data we do not hamper the bounds 
M\,M 00 ,M a . The most straightforward proof of local existence for smooth data yields a very 
short time of existence, but the a priori bounds permit the solution to be extended as long as 
(|125[) and (|128|) hold. Then we remove the regularization of the initial data, and pass to limit in 
the equation in distribution sense. We conclude that the limit obeys the equation, because the 
quadratic nonlinearity is strongly continuous in L 2 (for the right-hand side) and in W~ 1,2 for the 
advective derivative. Further details are left for the interested reader. 
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(130) G<\. 



We will pursue now the stability issue. Let us start by assuming that 
(129) —\\V x u{-Ml~^) < Ge- Kot 

Kq 

on the time interval t € [0,T]. We require that 
(130) 

Similarly, we assume 

(131) — [V x u(-,t)] a <He~ Kot 

K 

on the same time interval. We will show that we can choose G, H so that these assumptions are 
invariant in time: once initiated they cannot break down. Before we embark on the calculation, 
we should point out the simple strategy: we obtain from (|129l) . (|131[) a priori bounds on r. We 
then deduce a posteriori bounds for V x u that are strictly better than the assumptions (|1291 I131"j) . 
Because we work in a class in which the equation has unique solutions, and if the initial data for 
t are small, we deduce that, if initiated, the bounds continue indefinitely. 
We start by estimating the terms in (|66|) using (| 129[) and (|13ip . We will use 

(132) / 1 (z)dz < G 

Jo 

which follows from (|129[) . We have four terms, and in all of them, as well as in (|65l) we bound, 

/„' exp {~2 Ko (t -a) +2 J* j(z)} 7 (s)ds 

(133) < e 2G exp {-2n t} J* exp {2k s}k G exp {-K s}ds < Ge~ Kat e 2G 
< ^Ge~ Kot 

using (|132j) and the fact that AG < 1 which is (|130[) . From this calculation and (|65l) it follows that 

(134) \\T(;t)\\ L oo m < Ce' Kot [\\po\\ L °°^) + IItoIUoo^)] = Ce^M^ 

with C depending on d alone. Using (p2^1 and (I57|) in P5|) and integrating we obtain in the same 
way 

(135) \\r(;t)\\ LHmd) < Ce^ [\\p Q \\^m + Hh^)] = Ce'^M, 
We write the inequality (1551) as 

(136) [r(-,t)] a < I + II + III + IV 
and we start with 

(137) / < C[ Po } a e aG (I) e~ Kot = C[ Po ] a e- Kot . 



in which we used (|132[) and (|133[) above. This illustrates our use of the generic constant C: we 
used aG < \. For II we obtain 

(138) U < CWpoh-^e^e-^ J *[V x n(-, s)] a ds 

< C||p || L oo (R£i) e~ Kot J*[V x u{-, s)] a ds 

We used 2aG < 1 and (j!33[) . For /// we need to use the assumption (|131[) and instead of (|133[) we 
use a similar bound with the explicit KoGexp (— not) replacing j(t) and KoH exp (— ko<) replacing 
[\7 x u) a (t). We obtain: 

(139) III < C\\p \\ La , {Rd) e (2+2 ^ G He- Kot = C\\p \\ L ^ (Rd) He- Kot 
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Finally, 
(140) 



IV < c 



[ro]a + UrollioofRd) / [V x u{- , s)] a ds 



-2«ot 



where we use again (fl32l> and (2 + 2a)G < 1. Collecting terms we have proved that (fl29l [130} and 

([TBTjl imply that 

(141) 



[r(-,*)] a <Ce K °* |[po]e« + [To\a + #|| Poll (R^) + [ 1 1 PO 1 1 (K^ ) + ||n)|| L=°(R<*)] J [V x u(- , s)] a ds 

Now from (|87l) . and a Gronwall-like inequality we derive 
(142) [r(-,i)] Q < C [[pole* + Na +-Hl|po||z,°°(R<*)J exp<^ -t« + 



Indeed, (|14ip together with (|87l [93]) imply an inequality of the type 



y(t) < Ce~ tK ° 



a + b y(s)ds 



F(t) 



for the positive quantity y(t) = [r(-,t)] a , with 

a = [pola + [Tola + H\\p \\ L ^^d) 

and 

b = /c[||p ||L~(R d ) + ll T o||L~(Rd)]- 

Differentiating the right hand side F(t) and using the fact that y(t) < F(t) we deduce 

dF 



dt 

dividing by F and integrating we obtain 



log 



< -K F + Cbe- tKo F 



Ca 



Cb 

< -Kot H 

K 



and substituting back, we deduce (|142l) . Let us consider now the right-hand side of (EJT1) with 
(j = t. We recall that (fT3"I)) is 

(143) IK-.tJIUco^j^Ce-^Moo, 
(f!3"5|) is 

(144) HrMIUifR*) ^Ce-^Mi, 
and (|142p becomes 

(145) [r(.,t)] a < C[M Q + J ff||p ||L=o (R < 1) ] e - K0t exp{CDM oo } 
Assembling the right hand side of (|9Tj) we deduce 



(146) 



l^r| 



< CM^e"^' 11 + log 



1 | M* +a [Mp + gMpp] ' 

Moo 



Note that the term e K[)t cancelled in the expression for the argument of the logarithm. Now we 
use (|57f to obtain a posteriori estimates. From (|146p we obtain 

^||(VxU(-,i)|U°o(]R<i) 

(147) 



< CL»M oo e- K0 * U + lo, 



-, , M ° + ° [ Me, +HM oa ] 
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Also, from (|145|) and (|8"T|) we have 
1 



(148) 



KQ 



[(V x «)(-,*)]a < CD[M a + HMoo] exp{-K t + CDM^} 



H = 



Ma 



— [(V»u) (;t)] a < CDM a exp{-K t + CDMoo] 

Kq 



Let us choose G = 1/4 and 
(149) 

and so (|148[) becomes 
(150) 

The bound (I145[) becomes 

(151) [r(-, i)] Q < CM a e~ Kot exp {CM M } 
The upper bound (| 147[) implies the bound 

(152) — \\(V x u){;t)\\ La > m < B x Ge-^ 1 
where 

(153) B l = Bo 
and B Q = Bo (Mi, M a , M x , D) is given by 

(154) B (Mi, M a , Mqo, D) = CDM^ 



CD 2 Ml 



1 + log 



1 



M" 1 M Q d+c 



with appropriate constants C depending on d only. The functions Bo(M\, M a , Moo, D) and 
B 1 (Mi,M a ,M 00 ,D) are continuous in their arguments Mi > 0, M a > 0, M m > 0, D > 0, 
and hence locally bounded, and vanish when Moo — ► while the rest are held fixed, or when 
D — > and the rest are fixed. 

The bound (|150p implies the bound 

(155) — [(V s «) (;t)] a <B 2 He~ Kat 

Ko 

where the function B-i = B2(DM 00 ) is given by 

(156) B 2 = CDMoc exp {CDM^} 
with an appropriate constant C. We have proved therefore 

Proposition 3. Let r(x,t) be the solution of Jff^[ ) given in |^6j ) where p(x,t) solves V2b}) and 
assume that u obeys \129\) , ilSl]) ) on the time interval [0,T]. Then r and u obey the bounds \14^\ 
144\ M51\ \15'A \155\) one the same time interval, where the constants Mi, Mqo, M a are given in \96\ 



~97\ \98\) , the constants Bq, B\, Bi are given in \15tA \154\ \ 150(1 and the constants C are independent 
of solutions, fixed, and depend only on d. 

We combine now this proposition with the good local existence result, Theorem [T] We start 
with initial data that are small. To fix the notation, we denote by C± an upper bound for the 
constants C appearing in Proposition [3] above. We can arrange the constants so that, from 



(157) 
(158) 



[r(-,*)] a < TM Q , 
l|r(-,t)||Li(R4) <TMi, 
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and 

(159) ||r(^)|| L oo (Rd) KTMoo 
it follows that 

(160) — IIVxuO,*)!!^^) < — 

k v ' 8Ci 

and 

(161) -[^M;t)] a <^r 

because of the constitutive relations between u and r. Note that the size of r alone determines 
that of u. Indeed, in view of (1911). 



and in view of (193 



— \\V x u{-,t)\\L°°{R*) < CBq(TM\, TM a , TMoo, D) 
— [V x u(-,t)] a < CDTM a = CDYMooH 

So, the condition 

(162) B (rM x ,rM ,rM 00 ,.D)< 
implics (I160p and the condition 

(163) DrMoo < 

implies (|161|) . Then, we have automatically that (|129f) and ()131[) hold on the interval of time 
[0,T ]. In fact, if e KoT ° < 2Ci we deduce that 

(164) — ||V a u(-,t)|| L » (Rd) < iGe- Kot 
and 

(165) J_[ V:cU (.,t)] Q < ii/e-^ 4 

k 2 

hold on [0,T ]. On the other hand, as long as (fT29l[13T|) hold, the decay estimates (fl43 l HA4 \ [TgT| 
imply that 

(166) \\T(;t)\\ L ao m < |/l4o, 

(167) ||T(-,t)|| z i (Rd) < |m x 
and 

(168) [r(-,i)]«< ~M a 

provided e K °* > |Cie Cl£>M °°. The lower bound (|101j) on the life-span of the solution is 

(169) K Q T Q B Q {M x ,M a , M x , D) > e 
and it shows that if 



(170) B (Mi , M a , Moo , D) 



< e 
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then the solution exists for a period of time that is long enough so that r is reduced from its 
original size (|1661 11671 11681) at the end To of the interval. Requiring that 

'4 N 



(171) CxDM^ < log 



allows to arrange for T so that ^Cie ClDM °° < e K ° T ° < 1C\. Therefore, if the initial data satisfy 
(|1621 11631 11701 I171[) then the assumptions of Proposition [3] are satisfied with T — Tq. We can 
repeat now argument taking the same time step. Therefore we find that the solution, if unique, 
exists for all time and obeys the bounds (|157I1581 I159[) for all time. But then, inspecting the 
bounds (|65p and (|66j) in Proposition [T] we deduce that r decays exponentially, and consequently so 
does u. 

We note that the smallness conditions (|1621 11631 11701 II 71[) all follow from a single condition 

(172) B (M 1 ,M a ,M oo ,D)<e 1 

Theorem 2. Let the initial data r ,p satisfy (r ,p ) € (L 1 ^)) 2 n (C a {R d )) 2 , < a < 1. Let 
Mi, Mqo o,nd M a defined in Iff 61 \97\ \98\) denote the size of the initial data. There exists a constant 
Si such that, if 

(173) DMoo jl + log (l + M^M^MfA j < e x 

then there exists a unique global weak solution of (r,p) 6 L°°([0, oo), (L 1 (M d ) n C a (R d )) 2 ) n 
C([0,oo), (W~^ p (R d )) 2 ), p < oo of $2$, (9$, with u given by The meaning of weak solu- 

tions is in the sense of distributions, and the time derivatives of r and p belong to W~ 1,P (R ) for 
any p < oo. The gradient of velocity decays exponentially in time 

(174) — HV.uC.OIU-CB*) <Ke~* ot 

Kq 

1 



—not 



(175) — [V»u(-,t)] Q <Ke- K0 

The norms [r] a and || t|| decay exponentially: 

(176) [T(-,t)] a <Ke- K "\ 

(177) \\r(;t)\\ LP{Md) <K p e~^ 

for 1 < p < oo with explicit constants K that depend on d, D, Mi, Moo, M a , p and a > alone. 
The norms of p remain bounded. 

(178) l|p(-) i )llip(R d ) = IIPollif(Rd) 
and 

(179) [p(-,t)] a <K. 

Remark 1. The condition Ji7ff[ j is satisfied if the Deborah number is arbitrary, the initial tq is 
sufficiently small in C a (M. d ) and the initial pq is small in C a . Nontrivial initial data with small 
L°° norm and large spatial derivatives are allowed as well, because of the logarithmic dependence 
on M a . The conditions are also satisfied if the initial data are of order one but the Deborah number 
is small. Lf the initial data are smoother, the smoothness propagates. 

Remark 2. The proof of uniqueness, given below, uses Lagrangian transformations. Unlike the 
case of Euler equations, a proof in Eulerian coordinates seems difficult to obtain. 
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Proof. The previous argument is complete, modulo the uniqueness of the solution. We present 
here the proof of uniqueness. We consider the equation 

(180) jX = F[X] 

with X(a,0) = a, where X = X(a,t) is viewed as an element of X := C([0, T], C 1 ' a (M. d ) d ). We 
fix T. The function F depends on the whole path X not only on the value of X at some point 
and is obtained as follows. First we construct, using the initial data po G L 1 (R d ) n C a (M. d ), 
ooO) = Po(a)l + r (a) G (L 1 {R d ) n C Q (R d )) d2 , and using @2J|, the map X ^ a[X] given by 

rum 

1 j - 2 KoPo (a) Jj e-^-tya, i, a ) g (a, i, s) T ds + e - 2K °*(V a X(a, i))a (a)(V a X(a, i)) T 
where q(a,t,s) = q[X](a,t,s) = (V a X(a, t))(V a X(a, s)) _1 is given in p9]l. We take 

(182) T[X](a,t) = a[X](a,t) - ^(a)!, 
consider A(x,t) = X^ 1 (x,t), compose r[X] with it, 

(183) T X {x,t) =T[X](X- 1 (x,t),t), 
and solve the Stokes system, resulting in (|85t 

(184) u x (x,t) = khT 1 (RiTx + RiRmRnirT))- 
We write symbolically 

(185) u x = kA.- l m{T\X] o x- 1 ) = khr^rx) 

where EI stands for the combinations of Riesz transforms that appear in (|184[) 
(186) 

and thus fcA _1 lHI is the inverse of the Stokes system ((75)) . Finally, we compose with X(a,t) 

(187) F i [X](o,t)=u < x (X(o,t),t). 
Thus, is obtained via the succession of compositions 

(188) X i-> t[X] ^t x = t[X] o X- 1 hu x = kA^Wirx) i-> F = u x o X 
The norm in X is 

(189) \\X\\ X := sup \\X(;t)\\ C i,« m 

0<t<T 

We consider a fixed constant M and the set 

(190) V:={X eX | X(0,a) = a, - < det V„Jf(a, i) < ~, \\X\\ X < A/} 

The initial data for the PDE serve as parameters in the definition of F. We wish to show that two 
solutions X\ G T> and X2 G 2? of the equation (|180l) correspondingto the same po, tq, are identical. 
In order to do so we establish 

(191) \\{DF[X])Y\\ X <C\\Y\\ X 

with a uniform constant C that depends on M. We have to be careful to avoid taking derivatives 
of (To an d To- We start by noting that the map 

X h> a[X] 

is Frechet differentiable at X € T> as a map from the Banach space X to the Banach space 
£ = C(0,T; [C Q (R d ) n L 1 (]R d )] d ) of time-continuous maps with values in the space of matrices 
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with spatially Holder continuous and integrable coefficients. The derivative is a bounded linear 
map in C(X , E), 

Y i-> {Da[X])Y. 

The derivative has a complicated expression that depends on history but it can be easily obtained. 
The derivative of q[X] is 

(192) {{Dq[X])Y){a, t, s) = (d a Y(a, t))(V„X(o, s))" 1 - 9(0, t, s)(V a Y(a, s))(V a X(a, a))" 1 . 
Then the derivative of u[X] is 

((Da[X])(Y))(a, t) = 2 KQPo (a) /* e -2«o(t- s ) (i? g [jf]y)(o, i, s)g(o, t, s) T cfe 

(193) +2 Ko p («) /J e- 2K ^- s h(a, t, s){Dq[X]Y) T (a, t, s)ds 

+e- 2Kat [(V F(a, t))<7 (a)(V a X(a, t)) T + (V a X(a, t))a (a)(V a y(a, i)) T ] . 

It is clear that 

(194) |pa[X]F|| E < C||y|U 
where 

(195) IM|E= Sup {\\T(;t)\\ L ^ {Rd) + \\T(;t)\\ L i m + [T(;t)] a }. 

0<t<T 

The map 

X i-» t[X] 

is just a translation of <r[X] by an X-independent amount, so it has the same derivative 

(196) (Dt[X])Y = {Da[X])Y. 
Now, the map 

X(a,t) i-> T X (x,t) 

is obtained from the map t[X] by composition with X~ x . We cannot afford to take the derivative 
of tx- However, ux is obtained from tx using a linear smoothing operator of degree minus one, 
so we can differentiate it. We take a path X £ (a,t) = X(a,t) + eY(a,t), differentiate with respect 
to e and then set e = 0: 

(197) ((Du x )Y)(x, t) = kk- l W{{DT[X]Y) o X' 1 ) - kA- 1 B.((V x T X )(Y o X" 1 )) 

The first term appears when we differentiate t[X] o X^ 1 with respect to X , keeping X^ 1 fixed. 
The second term arises when we differentiate X -1 using the fact that 

±X-\ e=0 (x 7 t) = -(VxX-^YiX-^x^t). 

This fact is obtained by differentiating 4^ 

X- 1 (X(a,t) + eY(a,t)) = a, 

setting e = 0, and reading at a = X~ 1 (x,t). Then the derivative of t[X] o X^ 1 with respect to e 
is obtained using the chain rule 

£(r[X])(X-% =0 = ((V„r[I])oI-i)(|l-Vo 

= -((V a r[X]) o X-^iVxX- 1 )^ o X- 1 ) = -(V x t x )(Y o X- 1 ) 

The second term in (|197p is strictly speaking a distribution, as it involves derivatives of tx which 
is only Holder continuous. At this stage, we can view this as a formal calculation that will be 
justified in the end. Finally, we need to compose back with X(a,t). 

(198) (DF(X)Y)(a, t) = ((Du x )Y)(X(a, i),t) + {W x u x ){X{a, t),t)Y(a, t) 
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The second term in the expression above is unbounded as a linear operator from C 1,a to itself, 
simply because the coefficient V x u x is not differentiable. But this term combined with the second 
problematic term in (|f 97[) produces a commutator that is better behaved, AT[AT]y. Thus, we have 

(f 99) (DF(X)Y)(a, t) = k { A- 1 H((L>r[X]r) o A)} (X(a, t),t) + (K[X]Y)(a, t) 

where 

(200) (K[X]Y){a, t) = (V x u x )(X(a, t), t)Y(a, t) - k^U^V x t x )(Y o X- l )))(X{a, t),t) 

We would like to show that the map Y H> A"[X]F is a bounded linear operator in C(0, T; (C 1,a ) d ) 
with norm uniformly bounded for X G T>. Becuse composition with X^ 1 and composition with X 
are both bounded linear operators C 1 '" — > C 1,a , with norms controlled by M, the boundedness of 
K[X]Y is equivalent to the boundedness of the map 

(f)^>-L[X]4> 

where 

(201) <j>{x,t) =Y(X-\x,t),t) 
and 

(202) (L[X]4>)(x, t) = fcfVxA-^Tx t)4>(x, t) - fc(A- 1 H((V x r x )^))(a ; , t). 

It is important to specify the tensorial nature of this commutator, as not all such expressions are 
better behaved than their individual terms. In our case, Tx is a fixed symmetric matrix in E, k is 
a number, and the commutator is 

(203) {L[XW = k [cfPdvK-^mnTT - A^Hi^pdpTjT)] 
We can write the commutator as 

(204) L[X]<f> = k [^ p (d p A~ 1 M. imn )(rx n ) ~ (a p A- 1 H imn )(^rl")] - kK- x m imn {{d p ^)rT) 
For incompressible X we could only consider divergence-free <fr 

(205) d p f = 

but that would force us to work tangent to volume-preserving maps which would make the proof 
a little more complicated; we do not need to use that because the map 

J> ^ A- l W imn ((d p r)rD 

is bounded as a linear map from C(0, T; [C 1 ' a ] d ) to itself. Clearly, because tx is Holder continuous, 
and because classical Calderon-Zygmund operators are bounded in Holder spaces, there is no 
difficulty in bounding Holder norms of derivatives of the expression A~ 1 Mi mn ((d p (j) p )T x ln ) . Proving 
that the undifferentiated quantity is bounded is done using the fact that tx £ L 1 <~) L°°, and 
therefore (d p (jf)T X G L 1 n L°°. The operator A _1 H maps continuously L 1 n L°° to L°°. The 
operators TZ P i rnn = 9 p A _1 IIIi m „ are sums of classical Calderon-Zygmund operators and multiples 
of the identity. The commutators 

(206) ^ vn pmn {rT) - n pmn {^ P rT) 

are bounded as operators from C(0,T; [C 1 ' ]^) to itself. This is quite obvious for smooth tx but 
a little less obvious for tx e E. Let us write the kernel of lZ P i rnn as K p i mn , so the commutator is 

(207) / K{x-y){(j>{x)-^{y))T X {y)dy 

JR d 
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where we did not write all the indices and the time dependence for ease of notation. The kernel K 
is smooth away from the origin and is homogeneous of order —d. Differentiating in some direction 
and writing K' for the singular (of order d + 1) kernel obtained by differentiating K, we have 

(208) T[X\(j> = P.V.J K'{x - y){<f>{x) - (f>(y))rx(y)dy 

plus a nice operator (V x cf) p )lZpi rnn (T x ln ) ■ This last operator is clearly bounded in C a , with bound 
controlled by M, so we concentrate our attention on TLY]. Now we write 



(209) T[X]<f>= [ d\ 

Jo 

The kernel (x — y)K'(x — y) is homogeneous of order —d. It might have nonzero average on the 
unit sphere. Nevertheless, we subtract the value \7 x (f>(x): 



P.V. / {x-y)K'{x-y)V<t>(x + \{y-x))T X {y)dy 



(210) T 1 [X]4>= ( d\\p.V. [ 

JO I JR C 



(x - y)K'{x - y) {V x <j)(x + X((y - x))) - V x <j){x)} T X {y)dy 



The contributions left from the average on the unit sphere, if nonzero, are a constant multiple of 
(V x (f>(x))T X (x) and V x <f>(x)T 2 T X (x), 

(211) T 2 {t x )(x) = ( (x- y)K'{x - y){r x {x) ~ T X {y))dy 

both bounded with values in C a . The fact that Ti[X]<fr is bounded in C a , and similarly, that 
T2(t x ) is a Holder continuous function are classical. A proof can be found in ([3]). We have one 
more term in DF[X], namely 

(212) Y ^ kiA^MdDT^Y) o X- 1 )) o X 

Its boundedness is equivalent to the boundedness of the maps of the type 

(213) <j> i-> A -1 J!(ffxV(A) 

in C 1,a where g x is in S. These are easily bounded because when we take spatial derivatives we 
arrive at 

(214) h-> n{g x V4>) 

which are bounded in C a , and if we do not take derivatives, the L°° boundedness follows as above 
from the fact that g x V4> & L 1 <~) L°°. We have therefore verified the fact that DF[X] is bounded 
in X uniformly for X g T>. The function F is locally Lipschitz and because solutions of (|180l) start 
from the identity, they coincide for short time. The same argument does not need the initial data 
to be the identity, but rather the same invertible C 1 '" transformation so uniqueness propagates 
because F is locally Lipschitz. This concludes the proof of the theorem. 

5. A Regularization 

We consider here the system formed by the (JXJ) coupled with (|79p . (I5U1) via ([3]). In this section, 
we consider the same potential ([2]) but we allow R to be a function of x and t (but of course, 
not of m). The properties that a is symmetric, non- negative and the bound ((4]) remain valid. In 
order to see what dependences on R would be allowed by the energy considerations, we repeat the 
calculation leading to (|88|) : we take ([lj, multiply by fc(log /+{/), integrate, take (|79|) multiply by 
u, integrate, and add the two: we deduce 

(215) fc ^ ^ Rd ^ d f l0g fdmdx + 5 fvt d Tr ^ + k h$ T (/r< f Udm ) dx 
= - J Rd \V x u(x,t)\ 2 dx - ke J Rd J Rd f |V m (log/ + U)\ dmdx 
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It is thus clear that D t R > is energetically favorable. The solution on characteristics might be 
less explicit, however the equation obeyed by a is easily obtained by multiplying ([T]) by (m®m) / ' R 2 
and integrating. The result is very similar to (|29[) : 

26 26 ID R 

(216) D t a = (Vxu)ff + v(V x uf + ^ P l jj~ a 

We again see that D t R > has the effect of an additional damping. In fact, if is a constant, 
then the effect is precisely one of enhanced damping, and that is similar to the situation covered 
previously in Theorem [3J but in a better regime. More interesting, perhaps, is a damping that 
responds locally to very high rate of strain in the fluid. Let consider a coupled system in which, 
in addition to (P coupled to (|75|) , (1501) via © , R evolves according to 

(217) D t R = 5{\V x u(x,t)\)R 

and 5(g) is a smooth nonnegative function of one nonnegative variable g, that vanishes for 3 < § 

0, ifff<§ 



(218) 6(g) 



and satisfies 

(219) \S'(g)\<2C 
for all g > 0. In particular 5(|Va;u|) satisfies 



(220) 5 = Cq^/k 2 + \V x u\ 2 

if | V x u| > k > 0. The constant k > is fixed, and represents the order of magnitude of the largest 
permissible temporal growth rate. The constant Co > is chosen in function of the dimension of 
space (this is needed because of the tensorial nature of the calculation) so that 

(221) 3c\V x u\ - 2(5 < 

holds, if |Va;u| > k, where c is the norm equivalence constant that bounds max^ \diu l \ < c|Va;?x|. 

,2 Q 

Thus, c is a norm equivalence constant in R , and we can choose Co = §c. 

Because the problem evolves in time, the initial data /o and Rq need to be specified. We assume 
that /o is non-negative and smooth enough in x, decaying fast enough in m, (fo £ W 1,p (R d ; i 1 ((l + 
\m\ 2 )dm)) and R (a) > R mm > is in W 1 ' p (R d ) with p > d. Then, from pT6| and from the fact 
that a n are nonnegative for each fixed i, it follows that 

(222) D t Tr(a) < 2cn Tr(a) + 

R- 

which, after integration on characteristics, results in 

(223) supTr{cr(x,tj) < e 2 



supTr(CT (x)) H ^ — ||Po||z,<»(R<i) 



nun 



= N (t). 



In view of ((U), we thus have control of the L°° norm of a in time. Once this is achieved, based 
on previous results for similar models, we may expect to prove regularity. The method of proof of 
regularity given boundedness of a we employ here is the simplest and most explicit. The idea is to 
differentiate the equation, pay the price of differentiating the advective term, and do an Eulerian 
calculation in LP spaces, with p > d . Working in Lebesgue spaces makes it easy to take advantage 
of incompressibility. The approach used in the proof of global existence for Smoluchowski equations 
on compact manifolds coupled with time-dependent Stokes equations in ([6]), and in the proof of 
global existence ([11]) for Smoluchowski equations coupled to Navier-Stokes equations in d = 2, can 
be adapted to the present situation as well. In other words, we can consider the time depending 
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R evolving according to (|217[) . a Fokker-Planck equation ([T]) with potential given in ^ coupled 
via ([3]) to a fluid velocity evolving according to the time-dependent Stokes equation in d = 3 or 
Navier-Stokes equation in d = 2. The method of ([£.), used in (0) to prove the global regularity 
for d — 2 Smoluchowski equations on compact manifolds coupled with Navier-Stokes equations 
requires a little less smoothness on the Navier-Stokes initial data. That method is also Eulerian, 
but uses commutation and a penalty, that results in a controlled loss of regularity. The proof is a 
bit more technical as it employs Besov spaces and paradiffcrcntial calculus, but in principle it can 
be adapted to the noncompact particle phase space case with spatially depending R. We will not 
pursue these matters here, but rather content ourselves with the simplest proof, in the simplest 
nontrivial case. Differentiating (|2 16[) we obtain 



D t (d k *V) = -(9kU l ){dia») + (d lU l )(d k a^) + (d^)(d k a u ) - 25{d k u») 
(ZZi) +(flg I «*) £ r« + (0 fe W " (d k U)a^ - 2d k (6)(aV) + d^UpS.j) 



We mutiply by d k {a^) and sum. The terms involving explicitly first derivatives of u are bounded 
using (|22T1) : 

[-{d k u l ){dKjV) + (d lU %d k a^) + (d lU 1)(d k a il ) - 25(d k aV) - ^(d k a»)] (d k a^) 
< 3ck\V x (j\ 2 . 



(225) 



The terms involving explicitly second derivatives of u are bounded using (|223[) : 

(226) [(3u«V J + (d 2 k iu j )a il ] (fro**) < CN Q (t)\V x V x u\\V x <T\. 
The term containing a derivative of S is bounded using (|219l which implies 

(227) \V x S(\V x u\)\ < 2C \ 
and therefore 

(228) - 2d k (6)(a ij )(d k a ij ) < CN (t)\V x V x u\\V x a\. 



Summarizing, we have so far, pointwise: 
(229) 

-D t \y x a\ 2 < 3cn\V x a\ 2 + CN (t)\W 



■CN (t)e 



V 4 



IVxO-I +Ce 



X R 2 



Nx<r\ 



We have to bound the terms involving \7 x p and Vj-^j. In view of the fact that p(x, t) — po(A(x, t)), 
we have 



(230) 



R 2 



R 2 



(y a po)(A). 



Now, because of (1571) and (|217[) we have 

'V X A" 



R? 



(231) D t 
and, in view of (|221j) we deduce that 
(232) 

and consequently 
(233) 



R 2 



((V x «) + 251) 



\V x A{x,t)\ 
R{x,t) 2 



< 



Ri 



\VxP{x,t)\ <c e 



R(x,t) 2 



R 



2 

min 



\V aPo (A(x,t))\<C 



R 



2 

min 



|V a po| 
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The term involving W x (^2) is treated using (|217l) 
(234) Df ^ = (_ (V .„f-MI)(^) + 



and therefore, in view of ()221|) and (|227|) we deduce 



(235) 



D, 



\V X R\ 
R 3 



< CK 



R 3 



C- 



and, integrating on characteristics, we obtain 



(236) 



\V x R(x,t)\ 
R 3 (x,t) 



< Ce c 



\V a R (A(x,t))\ 
R 3 - 

mm 



r>2 



R 2 
R 2 

\V x V x u(X(A(x,i),s),s)\ds 



Now we collect the terms in (|229[) , divide by \V x a\ and use (|233p and (|236[) : 

D t \V x a(x,t)\ < 3cK\V x a(x,t)\ + CN (t)\V x V x u(x,t)\ + C^-e\V aPo (A(x, t))\ 
UM ' ! +Ce(N (t) + Wpoh-^e^ ( lY^^M + J* \V x V x u(X(A(x,t),s), S )\ds) 

The inequality f)23T|) has the form 



(238) 
where 

(239) 



D t {y{x, t)) < 3cKy(x, t) + C(t)z(x, t) + D(t) / z(X(A(x, t), s), s)ds + E(x, t) 



( y(x,t) = \V x a(x,t)\, 
z(x,t) = \V x V x u(x,t)\ 

C(t) = CN (t) = Ce 2 ^ [sup x Tr(a (x)) + -^-\\p \\ Laom 
D{t) = C^-(N (t) + || P0 || LO o (K£i) ) e c « t , 

^ E(x,t) = C^~e CKt {\S7 a po(A(x,t))\ + (N (t) + ||po|U~ (R< ^ l r '""' 



x,t))| j 



The interested reader can check that (|237j) is dimensionally balanced. Now it is time to start 
looking at LP norms. We start by noting that 



(240) 



J* z(X(A(x, t), s), s)ds 



LP(d x ) 



^ Jo P{X(A(x,t),s),s)\\ LP ^ds = f Q \\z(X(a,s),s)\\ LP{da) ds 
= Sq \\z(x,s)\\ LP{dx) ds 

because of incompresssibility. We integrate (|238[) on characteristics and take the LP norm. In 
order simplify the answer we use the fact that C(i), D{t) are non-decreasing functions of time. We 
obtain 
(241) 

\\y(-,t)\\ LP{Rd) < e 3cKt l\\yo\\ L p(M<i) + \\E(-,t)\\ LP{Rd) + (C(s) + D(t)(t - s))\\z(-, s)\\ LP{Rd) ds 
We recall ([57| . which implies 

(242) V.V^u = kR(V x a). 

In view of the well-known boundedness of Riesz transforms in LP spaces, we deduce that 

(243) ||z(-,t)IUp(R4) <C*||y(-,t)IUpCR««) 
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We note also that ||-E , ( - ,t)||LP(K< i ) i s explicitly a sum of norms of initial data multiplied by expo- 
nentials of time. Now a simple Gronwall argument provides an apriori bound for y in LP 

(244) h(;t)hp(M")<F P (t) 

with F p (t) an explicit function of time, with exponential growth, and depending only on norms of 
intitial data HooH^^), ||/3o|U°°(r<J), II V x a \\ lp (m*) , and Rmin, || V x -R |U*>(r<J). 

Theorem 3. Let f solve (QP with U given by (0) and R evolving according to ^217^ with smooth 
5 satisfying i219\) and H221\) . Let u be obtained by solving (79[) ; \80\) , with a given by (0). Assume 
that the initial distribution /o and Rq satisfy 



(245) 



sup xeR<) / Rd f (x, m)(l + \m\ 2 )dm < oo 
R (x) > R min > 0, 
/ Rd \V x R (x)\ p dx < oo, 
, Ik" (Uii 1 + \m\ 2 )\V x fo{x,m)\dm) p dx < 



with p > d. Then the solution exists for all time, is unique and obeys the a priori bounds 
and (Egp. 

Proof. The space W 1,p (M. d ) is a space of local existence and uniqueness of solutions. The condition 
m£ x R(x,t) > Rmin is invariant in time, because R(x,t) > Ro(A(x,t)) > R m in- Then the apriori 
bounds are enough to finish the proof. 

Remark 3. The a priori bounds hold in any LP , 1 < p < oo. 

Remark 4. Similar theorems hold if we replace steady-state Stokes equation with time- dependent 
Stokes equation in d — 2,3 and with Navier- Stokes equation in d = 2. 

Remark 5. Higher regularity of solutions can be obtained without difficulty from higher regularity 
of the initial data. 



6. Large data 

The problem of global regularity for arbitrary smooth initial data is open. The system formed 
by ([!]) and ([79)1 has potentially finite-time blowup. Indeed, consider the system in d = 1 



(246) 



D t f{x, m, t) + u x d m (mf(x, m, t)) = e<9 m (/<9 m (log f(x, m, t) + U{m))) 
U(m) = ^ 

a(x, t) = ^rf(x, m, t)dm 
u T = kHa 



where H is the Hilbert transform. This system is a 1-d analogue of ([TJ with @, (J2]) and (|57|) . 

2 

This blows up in finite time. Indeed, we multiply the linear Fokker Planck equation by ^ and 
integrate dm to obtain the analogue of ([2T)|) : 

(247) D t <j = 2u x a - + |j 
which then, in view of (|246|) is 

(248) D t a = 2koLLo - |ja + |j 
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and resembles the baby vorticity equation (l8j). The blow up argument of ([5]) works here also, 
notwithstanding the fact that the present equation is computed on characteristics. We form z = 
Ha + ia and use the fact that 

H(aHa) = \ {(Ha) 2 - a 2 ) 

to deduce 

, 2 2e .It 
DtZ = kZ ~& Z + l W 

which blows up in finite time. 

This simple example does not capture incompressibility and the tensorial nature of the problem, 
just like the baby vorticity equation of ([8]) does not capture them. The main quadratic nonlocal 
nonlinearity of (|29f is modeled by aHa and this is the only available model in one dimension. Unlike 
the incompressible Euler equation however, the system we study has an additional dissipative 
structure, and we believe that this is of some significance and represents the main fallacy of the 
one-dimensional model. We describe below a simple scalar model that addresses this issue. We 
take henceforth d = 2. We consider the variables 

( a(x,t) = \ (a n (x,t) -<r 22 (x,t)) , 
(249) I b{x,t) = a 12 (x 7 t) = a 21 (x : t) 7 

{ c(x,t) = a n (x,t) + a 22 (x,t) = Tr(a(x,t)) 

The gradient of velocity is represented by X(x,t), n(x,t) and uj(x,t) given in (fS§|) . (|70j) and (jTTj) . 
The equations (|2"9"|) can be written as the system 



(250) 



D t a = —ujb + Ac — -j^a 
D t b = uia + /ic — -j^b, 
D t c = 4Aa + 4/x6 - j^c + A-^p 
D tP - 



As we saw before, if we couple this with an equation for the velocity (steady or unsteady Stokes, or 
Navier-Stokes) , then the regularity issue is decided by whether or not we can bound c in L°° (dx) . 
We note in passing that the co-rotational system corresponds to A = /i = in the system above, 
and the bound for c follows immediately. Let us multiply the c equation by |, the a equation by 
2a, the b equation by 2b and subtract the last two from the first. We obtain 

c 2 2 , \ 4e f c 2 , , \ 2e 



(251) D ^__^_^ = __^__^_^ + _ pc 

This cancellation of nonlinearity is not surprising because 

c 2 

(252) — - a 2 -6 2 = Det(a) 

and the determinant is conserved along particle trajectories if e = 0. The matrix a is symmetric 
and positive by construction, and is given in terms of a, b, c by 

(253) '=( f r f-a 
The positivity of the matrix is equivalent (in this case) to the positivity of the determinant, i.e. to 

c 2 

(254) — - a 2 - b 2 > 0. 
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The two eigenvalues of cr 

(255) «i, 2 = °- ± ^a 2 + b 2 

are both positive. We have of course c = z\ + z 2 > and z\ — Z2 = 2V<x 2 + b 2 - Because of (|254[) . c 
controls Vet 2 + b 2 . But, on the other hand, growth without bound of c on any trajectory, cannot 
happen without growth without bound of \/a 2 + b 2 on the same trajectory. Indeed, if a particle 
path would be such that \J a 2 + b 2 is bounded on it, but c grows without bound or blows up in 
finite time, then, for large enough time we would have 



on the particle path (because p is bounded) and then, from (1251 1) we would arrive at the contra- 
diction that c remains bounded. From (|250[) we can write 

(256) c(x,t) = e-&*co(A(x,t)) + J e -^ (t ~ s) (Xa + fib - -^p ) (X(A(x,t), s), s)ds 

Thus, exponential growth or blowup of c is possible only if (Xa + pb) has time integrals on particle 
paths that are positive and grow exponentially or stronger, without bound. 

In two dimensions we can express the velocity in terms of a stream function ip(x,t) and write 
it 1 = —d2ip, u 2 = diip. Then uj — Aip, and therefore 

(257) A = ft.Sb(-A) _1 w = Bu 
and 

(258) fi = -^(d 2 1 -d 2 2 )(-A)- 1 u J = -Au J 
where 

(259) A= l -(dl~dl)(-A)-^\{R\-Rl) 
and 

(260) B = did 2 {-A)- 1 = RxR 2 

are bounded operators in L P (M. 2 ) spaces. They are also bounded, selfadjoint in L 2 (R 2 ), they 
commute AB = BA, and each is given by a multiplier at Fourier level, 

(261) Th(0 = 4^r^) 

and 

(262) Bh(0 = -^-h(0. 
Note that 

(263) A{A 2 + B 2 )=1 

Let us consider now the time independent Stokes system ([79]) . Taking the curl of (|79"1) and inverting 
the Laplacian we obtain 

(264) uj = 2k (Ab - Ba) 
Consequently, from (|257|) 

(265) A = 2k (-B 2 a + ABb) 
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and from (|258l) 

(266) fx = 2k (ABa - A 2 b) 

It is convenient to measure time in units of wr and then the system (|250|) coupled with (|264l) , (|257l) 
and (|2"55|) is 

{D t a = -(A(b) - B(a))b + [-B 2 (a) + AB(b)]c - e~ x a 
D t b = (A(b) - B(a))a + [AB(a) - A 2 {b)]c - e~ x b 
D t c = A{[-B 2 {a) + AB{b)]a + [AB{a) - A 2 (b)]b} - e~ l c + 2e~ x p 
D tP = 

where e is given by (|99l) . The all- important term Xa + ph is given by 

(268) — (Ao + fib) = aB(uj) - bA(u) = -aB{B(a) - A(b)) + bA(B(a) - A(b)) 

This expression is quadratic, nonlocal, and has negative spatial integral. While blow up requires 
the pointwise positivity of this expression (|256[) . its spatial average is negative. Integrating the 
third equation in (I267P we obtain, using the selfadjointness of A and B: 

(269) Jr 2 c (f ' V dx + 4 /o \B(a)(x, s) - A{b){x, s)\ 2 dx 
= / R 2 c (x)dx + 2(1 - e~") J R2 p (x)dx 

This is just the energetic bound on J Tr adx. We saw earlier that, if the size of initial data is of 
order one, then small enough e leads to global existence of solutions. The case of large initial data, 
moderate e is wide open. In order to clarify the issues, we will take e = oo. In this case the system 
is simplified somewhat, because, in view of (|251l) we have 



(270) c = 2^ a 2 + b 2 + d (x,t) 
where 

(271) d (x,t) = Det(a (A(x,t))) 

and so, the first two equations of (|250|) become just 

D t a = ~b{Ab - Ba) + 2^(a 2 +b 2 + d ) B(Ab - Ba) 
D t b = a(Ab - Ba) - 2y/(a 2 + b 2 + d ) A(Ab - Ba). 



(272) 
and consequently 

(273) D t y/(a 2 +b 2 + do) = 2aB(Ab - Ba) - 2bA(Ab - Ba) 

Integrating in space we again get the bound on the trace. Blowup for the system (I272p is still too 
difficult to analyse. The simplest possible didactic model of a scalar equation exhibiting this kind 
quadratic nonlocal structure with an L 1 dissipation is a scalar equation of the form 



(274) d t r = -t{A 2 t) 

where A is a bounded selfadjoint operator given in Fourier representation by multiplication by a 
function that is homogeneous of degree zero, like above, 

(275) Th(0 = a(0H0 

with a(A£) — a(£), a(£) € R, L 2 £*(£)d£ = 0, a smooth on the unit sphere. The unknown r(x,t) 
(representing the trace of a) is a positive scalar. We note the salient features of this model. If tq is 
smooth and positive, (let us consider for instance r(-, 0) € L 1 n C°' s with s > 0) then the solution 
exists and is unique for some time, and remains positive as long as it exists. If a bound on ||t||l«. 
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is given, then the higher regularity of the solution follows. Integrating the equation, we have the 
dissipative bound 

(276) J r{x,t)dx + \\AT\\ 2 L2 dt< J T{x,0)dx 
and multiplying by A 2 t and integrating we have 

(277) / \Ar(x,t)\ 2 dx < [ \At(x, 0)\ 2 dx 



We are going to give an example of global regularity for some large data for such an equation, in 
the spatially periodic case. For simplicity we take the period to be 27r in each direction. We write 

(278) T (x,t)=T (t)+ J2 T k(t)e ik - X 

fcez 2 \{o} 

The requirement that t be real is implemented by t% — T_fc. We will consider the symbol a{k), 
and assume that it is real valued, even, a(— k) = a(fc), and bounded |a(fc) | < T. We also assume 
a(0) = 0. The equation (|274|) is the infinite system of ODEs 

(279) ^ = - E 

k+j=i 

Proposition 4. Let t(x, 0) = J^k T~k(0)e x satisfy 



(280) r_ fc (0) =r fc (0), 

(281) ^(1 + W) S M0)| <C S (0) <oo 
for some s > and 

(282) r (0)>^|r fe (0)|. 
Then the solution of \219^ exists for all time, and obeys 



(283) r_ fc (t) = r k {t) 

(284) ^(l + \k\Y\T k (t)\ < C s (0)e r+lTai0)T2t < oo 

k 

and 

(285) r (t) > \r k (t)\ 

fe#0 

which of course implies that t(x, t) remains smooth and positive. 
Proof. We start with the dissipation, the equation at I = 0, which reads 

(286) ^ = -E«(*) a l'*(*)l a 

k^O 

Now for I we have 

^jT = - X! T k(X 2 (j)Tj - T a 2 {l)Ti 

k+j=l, fc^O, j#0 
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and therefore 

(287) M <- ro a 2 (0h|+ £ M^CjJM 

k+j=l, kjtO, 

Summing in / 7^ we obtain 

iE M o hi < -ToE^ " 2 (Ohl + E^oE* +i =i, mo,^o W« a (?')N 

= - r oEi^o a 2 (0M +E^oE^o a2 0') kill 7 *! - Ej#o a2 (i)l r -ill r il 

= (E^rXWN) (-to + E^o W) + 3h> 

This results in the inequality 

(288) |^ h |_ ro j < ^^OOlr.l] [^h|-ro 

If the initial data is non-positive, the quantity remains non positive. Thus we have the invariance 
of this cone in function space. Once this is achieved, we know 

(289) ^h(t)|<T (t)<T (0) 

which implies an L°° bound on r. We know this should be sufficient for regularity. In the present 
situation, the proof of regularity is quite straightforward: we multiply ()287[) by (1 + |Z|) S and sum. 
Using |/| < 2max(|fe|, |j|) and \a(J)\ < T we obtain 

^E^ + i^n^ 2 ^ 2 (En 1 lE^+i^N 

Using (|289[) we deduce 

(290) |^ (1 + | Z | rN < 2 ^ ro (o)r 2 

which finishes the proof. 
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